Does the crossover from perturbative to nonperturbative physics in QCD become a 

phase transition at infinite N ? 
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We present numerical evidence that, in the planar limit, four dimensional Euclidean Yang-Mills 
theory undergoes a phase transition on a finite symmetrical four-torus when the length of the sides 
/ decreases to a critical value l^. For I > Ic continuum reduction holds so that at leading order in A'', 
there are no finite size effects in Wilson and Polyakov loops. This produces the exciting possibility 
of solving numerically for the meson sector of planar QCD at a cost substantially smaller than that 
of quenched 517(3). 



1. Introduction. Nonabelian gauge theories in 
four dimensions interact strongly at large distances and 
weakly at short distances. The major conceptual achieve- 
ment of lattice gauge theory has been to show that the 
continuum limit contains both regimes and interpolates 
between them in a smooth manner. At present, only nu- 
merical simulations can bridge these two regimes. 

It has been a long held hope that the task would 
simplify at infinite number of colors N. In its most opti- 
mistic version, the hope is that somebody will come up 
with an analytic solution to the N = oo, "planar", limit 
at all scales, and that it also will be feasible to compute 
-i- corrections. We are far from attaining this goal, but 
there has been recent progress on similar, albeit vastly 
more constrained theories the most notable example 
being TV = 4 SUSY YM. The exact exphcit solution of 
the planar limit in this case is limited to extreme cases, 
the simplest among them being when the 't Hooft cou- 
pling (A = 9ym^) is taken to infinity. The general case 
has been mapped into a two dimensional field theoretic 
problem which is not yet fully understood. The correc- 
tions will come from the interactions in IIB string theory 
expanded around an AdS^ x 5*5 background stabilized by 
a non-trivial RR flux. At large A many answers can be 
obtained by relatively simple calculations of small per- 
turbations around this supergravity background. 

In pure YM, we do not have a free coupling constant: 
instead there is nontrivial scale dependence reflecting the 
breaking of conformal invariance. A simple way to visual- 
ize the situation is to consider Bjorken's "femto-universe" 
[3 , where one studies the Hamiltonian of QCD restricted 
to a small three dimensional box of side I. Equivalently, 



one can not only shrink the system (which takes care of 
linear momentum) but also raise its temperature (which 
deals with frequencies). The latter situation is best de- 
scribed by Euclidean field theory on a four torus. A rem- 
nant of Lorentz invariance is preserved when the temper- 
ature and linear size of the torus are related so that in Eu- 
clidean space the four torus has sides of equal length. As 
I is varied, the bulk physics of the system is perturbative 
for IAqcb << 1 and nonperturbative for ZAqcd >> 1- 
Thus, the role of the coupling is taken up by L 

In this paper, we consider SU{N) pure gauge theory 
in the planar limit on a four-torus of side I. The most 
basic question is whether, as a function of I, this sys- 
tem undergoes any phase transitions at iV = oo. There 
can be such transitions because N equals infinity: for 
finite N all these transitions will be smoothed out into 
crossovers. It is possible that when the torus is taken 
to infinite four-volume in a specific way, the crossovers 
become transitions even at finite N. How the various 
transitions merge into a coherent picture is a question 
that we have only begun to explore. 

Recent work in three dimensions Q has shown that 
when / is decreased from infinity, there will be a transi- 
tion at I = Ic- For I > Ic, the system realizes a continuum 
version of Eguchi-Kawai reduction , whose salient fea- 
ture is that expectation values of arbitrary Wilson loops 
are exactly /-independent once N — oo. This continuum 
Eguchi-Kawai reduction is a property that is natural for 
a system of free strings but more difficult to understand 
in field theoretical terms where, in order to erase the 
perturbative /-dependence, one is required to enlist an 
averaging over a moduli space of minima of the classical 
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action. In this paper, we show that essentially the same 
effects work in four dimensions as well. 

2. General phase structure. We used a single pla- 
quettc Wilson action as the cost of simulation increases 
with the number of colors N as N^, making the sim- 
plicity of the action a relevant resource consideration. 
We worked on tori whose sides consist of L lattice sites. 
The physical side length is / = aL, where a is the lat- 
tice spacing. L was varied from 1 to 10. The parameter 
N was chosen to be a prime number, taking the values 

= 23,31,37. The preference for prime N is because 
the transition we found has to do with Z{N) groups, and 
for prime N, Z{N) does not have subgroups that could 
confuse the picture. In one case, we did a simulation at 
N = 27 and found nothing unusual, so it is possible that 
the choice of prime A'' values was unnecessary. 

E Tr[U,A^) + UlAx)] (1) 

U^A^) = U^{x)U,{x + ^)Ul{x + v)Ul{x) (2) 

We define b = jj = and take the large N limit with 
h held fixed. As usual, h determines the lattice spacing 
a. The lattice is a symmetric torus of side L. The gauge 
fields are periodic, x is a four component integer vector 
labeling the site, and fx either labels a direction or denotes 
a unit vector in the /i direction. The link matrices f/^ (x) 
are in SU{N). 

There is a Z'^{N) symmetry under which 

U^{x)^e^U^,{x) (3) 

for all X with a:^ = c^. The integers are fixed, and 
the integers label the elements of the group; c^. A;,, = 
0, 1, ..,L — 1. Changing the c^'s amounts to a local gauge 
transformation. 

Polyakov loops are denoted by -P|i(x) and defined by: 

P^{x) = U^{x)U^{x+fi)U^{x+2tj.)..U^{x+{L-l)fx) (4) 

Under the above symmetry, gets multiplied by a 

phase. The gauge invariant content of P^{x) is its set 
of eigenvalues (the spectrum) e*^'^, i = 1,2..., N. The 
ordering is not gauge invariant, and there is a constraint 
that detP^(a;) = 1. Under a Z{N) transformation, the 
set of eigenvalues is circularly shifted by a fixed amount. 
The spectra of P^ix) and of P^ix + jfj) are the same 
for all j = 0,1,2...., L — 1. Wilson loops are defined 
similarly to Polyakov loops, only they are invariant under 
the Z'^{N). Often we shall speak about the angles Of, 
thinking about them round a circle and referring to them 
also as the "spectrum" . 

At a given L, we increase b gradually, until a point is 
reached where one of the four Z{N) factors, acting in a 
randomly picked direction n, breaks spontaneously. The 



breaking is reflected by a change in the spectra of Pf^ (x) 
away from a form symmetric under circular shifts. This 
happens when a gap larger than ^ opens up in the the 
angle spectrum at some random location round the circle. 
This event can be detected in various ways. 

At infinite N, six phases are encountered as b is var- 
ied from zero to infinity on a lattice of size so long as 
L > 9. In the range < 6 < 0.36 the system is in a "hot" 
phase (denoted by "Oh"), where the Z^{N) is preserved 
and the 1x1 Wilson loop has no gap in its spectrum. As 
b increases one goes into a "Oc" phase, where the Z^{N) 
symmetry still is preserved, but the 1x1 Wilson loop 
now has a gap in its spectrum. The gap is centered at the 
point —1 on the unit circle, so that charge conjugation is 
also preserved. Next one goes into a "Ic" phase, where 
exactly one factor of the global Z'^{N) is broken. As b in- 
creases further, additional factors of Z'^{N) successively 
break until the phase "4c" is reached, which extends all 
the way to 5 = oo. For 5 < X < 8 the phases "Oc" and 
"Ic", including the transition between them, can be ex- 
tended downwards in b, as metastable phases, into the 
"Oh" phase region. Thus, using metastability, we can ex- 
tend most of the phase structure of interest from L > 9 
to L > 5. For L < 4 the stable "Oc" phase is "squeezed" 
out, and we are left with only five phases. The case L = 1 
is the original Eguchi-Kawai model. All this holds with 
Wilson's single plaquette action and might change with 
a different lattice action. However, physical properties 
that survive the continuum limit should be insensitive to 
the choice of action. 

The "Oc" phase is the most interesting phase because 
there planar QCD exhibits confinement and stringy be- 
havior at large distances and field theoretic asymptotic 
freedom at short distances. The phase "4c" is the phase 
where planar QCD is well described by Bjorken's femto- 
universe heated to high temperature. 

There is little doubt that the "4c" phase survives in 
the continuum limit. This means that there exists a finite 
range of torus sides I between zero and some small scale 
where planar continuum QCD is in a "4c" phase. There 
also is little doubt that the "Oh" phase does not have a 
continuum limit, i.e. there is no finite range of ^-values in 
which the continuum system is in a "Oh" phase. In otlicT 
words, the "Oh" phase is a lattice artifact. In this paper, 
we shall present evidence that the "Oc" phase does have 
a continuum limit, describing the system in the range 
oo > I > Ic. More work will be needed to complete 
the continuum phase diagram and see in detail how the 
system goes from the "Oc" phase all the way to the "4c" 
phase as the torus is shrunk in size. 

3. Numerical method and results. We sim- 
ulated the system using the Monte Carlo method em- 
ploying heat bath updates and overrelaxation updates. 
The heat-bath updates amounted to sequential SU{2) 
updates, going over a set of of SU{2) subgroups 

identified by choosing two distinct integers between 1 and 
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N |g. For most of the values of b we used, each SU{2) 
update was done using the Kennedy-Pendelton method 
|8j . For few smaU values of b we used the original Creutz 
method 0- The cost of a heat-bath update goes as 
as N increases. 

The overrelaxation update was a full SU (N) update 
Hol and had a comparable cost. Our implementation 
went as follows: 

The portion of the action S that depends on a partic- 
ular link matrix, denoted by U, Sr(U), is given by 



Tr [f/E] = ^[Sr{U)+iSi{U)] 



(5) 



with real Suj{U). S are the "staples" , a positive number 
(the coupling) times a sum of simple unitary matrices. S 
is determined by U and when this is not evident from the 
context we shall use the notation !][/. With probability 
one, S has non-zero determinant permitting a unique def- 
inition of a unitary matrix Vs; 
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S, detV^E = 



detE 
IdetEl 



(6) 



Vs is calculated as follows: EE^^ and E^E are diago- 
nalized using the Householder method. The eigenvalues 
are distinct with probability one and make up a positive 
diagonal matrix D. The diagonalizing matrices provide 
two representations of D: D = XT,'<T,X^ = FEE^rt 
and, finally, we end up with Vy: = X . 

Taking $5] to obey tt > $5] > — vr we define a new 
SU{N) matrix by 



V 



(7) 



The update starts by "offering" the replacement U 
V. The new action is given by 

SniV) = cos (^) Sn{U) + sin (^] Sj{U) (8) 



V N 



N 



For large TV we expect Sr{U) to be order N'^ and Si{U) 
to be order N. This implies that the change in action is 
order one and hence there is an order one probability 
of making a large move in configuration space. In the 
Metropolis step, the a priori probability for change is 
unity for the U ^ V transition and zero for anything 
else. Applied twice, this transition becomes the identity; 
therefore the Metropolis step only depends on R, the 
ratio of the Boltzmann factors. 



R — e 



(2ia)s,(£/)-[i-co.(2£a)]s«(C/) 



(9) 



The acceptance probability for the change is taken as 
min{l,i?}; this satisfies detailed balance. 

We found that the acceptance rate for the overrelaxed 
update was over 95% for all our N and b values. We em- 
ployed a mixture of heat bath and overrelaxation steps 



L=6 N=31 b=0.351 started from 1c 

(measurements taken every 50 lattice passes) 
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FIG. 1: History of the variable p{P^) for each direction. We 
see the evolution from a state where one of the four Z{N) 
factors is broken to one in which all four are preserved. 



of equal amounts. A more comprehensive independent 
study of full SU (N) overrelaxation has been recently pre- 
sented in 

Our Polyakov loops (as well as various Wilson loops we 
looked at) were built out of U^{x) matrices, rather than 
the original link matrices U^{x). The U^{x) matrices are 
defined in term of the U^{x) by an iterative "smearing" 
procedure jT3 | . One step in the iteration takes one from a 
set ujr\x) to a set uj:i^^~^^\x) , by the following equation: 

4"+i)(2;)^aC/(")(x) + i^E,,,„„ 



6 u'^^'ix) 



C/("+i)(:r)=4"+i)(x) 



1 



^[x^-^'\x)]^xir^'\x) 



We chose a = 0.45 and iterated L-times: 



U,ix)^Ul,^\x) 



(10) 



(11) 



This has little effect on Polyakov spectra at smaller b 
values, but, after the transition the well known ultravio- 
let renormalization j3| of Polyakov loops will reduce all 
traces TrP^(x, /z), effectively making the angle spectrum 
look more uniform and making the transition harder to 
discern. This effect is reduced by the above smearing. 
In our three dimensional work |j| we could do without 
smearing, since the ultraviolet divergence is milder. The 
Polyakov loops in terms of the smeared links have the 
same symmetry properties as the Polyakov loops in terms 
of the original links and therefore provide perfectly ade- 
quate order parameters for Z{N) symmetry breaking. 

We looked at several observables. Two were the most 
useful for identifying the phase transition. The first is 
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L=7 N=31 b=0.3568 started from Oc 

(measurements taken every 50 lattice passes) 
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FIG. 2: History of the variable p(P^i) for each direction. We 
see the evolution from a state where all four Z{N) factors 
are preserved to one where one factor is broken. During the 
first fifty passes (before the first measurement) Polyakov loops 
in direction 3 have acquired some structure but, ultimately, 
direction 2 is selected for breakdown and the Polyakov loops 
in the other three directions converge to a symmetric state. 



L=9 N=31 b=0.3660: one broken Z(N) 

Unnormalized distribution of maximal angle spacing 
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FIG. 3: Here we show the difference between the distributions 
of the largest inter-angle spacing for smeared Polyakov loops 
in different directions in the phase where exactly one Z{N) 
factor is broken. (At other couplings, where no Z{N) factor 
is broken, all four distributions look like the three unbroken 
ones here). 



FIG. 4: Angle distributions in four directions in the Ic phase. 
There are twenty seven periods in the superposed oscillations. 
The peaks, except close to the gap associated with direction 
3, are equally spaced. 

4D 2-loop P-function for L^(b) 



Tadpole Improved 



L„'=0.275 

L„'=0.260 

L„'=0.245 

♦Oc <-> 1c (stable) 

■ Oc <-> 1c (metastable) 

■ lc(metastable) -> Oh(stable) 
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FIG. 5: The transition ranges compared to possible two loop 
renormalization group curves with tadpole improvement. 

PiP^) = ^(Y.'^^ 2^^^-€f) (12) 

The averaging is over the 3-plane perpendicular to /i and 
over configurations. Equally spaced angles respect the 
Z{N) symmetry in the /i direction and maximize p to 
0.5. When the angle-spectrum starts getting modulated 
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TABLE I: 


Summary of ranj 


jes for the "Oc" to "Ic" transitions. 


L 




(^min ; ^max) 


5 


31 


(0.3470,0.3480) 


5 


41 


(0.3473,0.3485) 


6 


31 


(0.3510,0.3520) 


7 


31 


(0.3560,0.3568) 


8 


23 


(0.3590,0.3610) 


8 


31 


(0.3595,0.3605) 


9 


23 


(0.3630,0.3655) 


9 


27 


(0.3635,0.3650) 


9 


31 


(0.3630,0.3660) 


10 


23 


(0.3662,0.3678) 



and opens a gap, p drops below 0.5. The second observ- 
able that we found useful, p(sniax), is constructed as fol- 
lows: Among all spacing between adjacent angles round 
the circle select the largest one, Smax- Its probability 
distribution, /^(sniax), strongly depends on whether the 
Z{N) associated with the direction under consideration 
is broken or not. When the Z{N) is not broken the distri- 
bution just reflects universal angle repulsion. However, 
when a gap opens in the spectrum, it dominates the dis- 
tribution. Of course, looking directly at the histograms of 
the angles associated with individual directions remains 
the most direct way to observe the behavior of the sys- 
tem. Figures ^ and [21 show examples of the evolution 
of p when the run started from a typical configuration 
in the wrong phase. Figure |21 shows an example of the 
difference between the distribution of the maximal level 
spacing in the direction corresponding to a broken Z{N) 
and the maximal level spacing distributions in the other 
directions, whose corresponding Z{N)'s are unbroken. 

Figure 01 shows an example of the angle distributions 
associated with the four independent directions when the 
system is in the "Ic" phase. The superposed oscillations 
reflect the global SU{N) constraint on angle values. For 
a SU{N) matrix drawn with Haar probability measure a 
relatively straightforward calculation gives, with ^ — —'■ 



(-1) 



N 



■ cos(A^7r^) 



(13) 



Thus, the expected swing between minima and maxima 
should be ^ in our plot and this fits more or less. Also, 
because 27 is an odd number, we have a maximum, rather 
than a minimum, at ^ = 0. 

The "Oc" to "Ic" transition also breaks hypercubic in- 
variance since one direction is randomly selected by the 
breaking. The "Oc" to "Ic" is transition most likely is of 
first order, and therefore generates hysteresis cycles for 
relatively short runs. Our runs were of the order of few 
thousands and did not allow a very precise identification 
of the location of the transition or a definitive determi- 
nation of its order. We ran hysteresis cycles looking for 



two extremes which determine the range we believe the 
true transition is in. At the first extreme we start from 
a configuration typical of the "Ic" phase, and see that 
the system disorders, restoring the remaining Z(N) fac- 
tor. This is what happens in Figure ^ For our larger 
volumes {L = 9 and L — 10) we try to find the largest 
b where the configuration evolves in this way within at 
most 3,000-4,000 passes over the lattice. At the other 
extreme we start from a totally symmetric configuration 
and observe the system evolving into a "Ic" phase, like 
in Figure 121 Here we try to find the smallest b where this 
scenario is realized. The ranges in b so obtained were of 
length somewhere between 1 x 10"'^ and 3 x 10""^. A few 
checks, performed by varying N at fixed L , showed that 
the finite effects were at most of order 5 x lO^''. With 
this accuracy, we were able to check whether the location 
of the transition bc{L) varies with L in a way compatible 
with asymptotic freedom. Our numerical work used up 
about one year's worth of time on a dedicated desktop 
PC with a modern processor and 2GB of memory. With 
this rough map of the phases in place, one could proceed 
to finer determinations, but this would require one or two 
orders of magnitude more computer time. 

If the transition we are searching for truly is a con- 
tinuum phenomenon, the inverse of the function bc{L), 
Lc{b), should behave for & ^ oo as: 



Lc{b) 



11 



(14) 



The asymptotic regime is not reached at L ~ 10, but by 
the "tadpole" replacement 



bj = beib) 



Lcib) 



<b) = ^{TrU^^^ix)) (15) 



11 



121 2i,T^bj 

e 11 



(16) 



the asymptotic behavior is supposed to set in much ear- 
lier hj. The numerical effect the replacement of b 
by 6/ has is summarized by the approximate relation 
Sbj ~ 1.3(56 which holds in the vicinity of the transition 
at L = 9. With b replaced by 6/, the theoretical curve 
becomes somewhat steeper. The plaquette expectation 
value, e(6), is taken on the symmetric side of the tran- 
sition and using the MC data, can be well fitted in the 
range of interest by: 



e(5) 



1 



b 



+ 



1 I a2 I 03 



(17) 



When 6 ^ oo, e(6) = 1 at leading order in j and bj ~ b. 



■1 ana oj 

We varied L between 4 and 10 and b between .344 and 
.366. e{b) was reasonably well fit by ao = —0.58964, oi = 
0.08467, a2 = -0.50227, as = 0.05479. Hence, 



Lo 



0.1524L^ 



(18) 
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The ranges we determined for the "Oc" to "Ic" tran- 
sitions are reasonably well described by a range of Lq 
constants between 0.245 and 0.275. Figure [3 shows the 
ranges we established on a plot together with lines rep- 
resenting the tadpole improved two loop renormalization 
formula with different amplitudes. The relative consis- 
tency of this fit constitutes our numerical evidence that 
the transition is physical, occurring in the continuum at 
a finite scale. The relevant numbers for the "Oc" to "Ic" 
transition ranges that went into figure El are collected in 
Table D 

Let us first discuss all the stable phases, ignoring the 
metastable ones. In addition to the "Oc" to "Ic" transi- 
tion we have been focusing on, the system also undergoes 
a lattice transition from "Oh" to an "Xc" . For small vol- 
umes, X will be large than 0, but, starting with L = 9, 
X= 0. The "Oh" to "Xc" (X^ 0) transitions at smaller 
L values are strongly first order and occur at b = 6bulk- 
A similar transition also occurs at L = oo for any N > 5 
[l^ . In this case this bulk transition is well known to be 
associated with a large jump in e{b) at &bulk- At finite 
N the transition does not break any symmetry. 

At = oo, and any finite but large enough L, the loca- 
tion of the bulk transition can be estimated with the help 
of to occur at 6bulk = 0.3600. (For an earher deter- 
mination, see 01). So long as L > 8, the bulk transition 
again breaks no symmetry: It takes the system from the 
"Oh" phase to the "Oc" phase. Again, e{b) undergoes a 
significant jump at &bulk- In addition, with TV = oo, the 
average eigenvalue distribution of the plaquette variable 
now undergoes also a qualitative change, opening a gap 
at angle tt when b increases through ^bulk- 

We have investigated the L = 1 case in some detail. 
This case is special, and the algorithms we used are some- 
what different; since this is a bit of a side issue, we shall 
not elaborate in detail. We found evidence for five stable 
phases. Examples of firmly determined couplings in each 
phase are: 6=0.150 in "Oh", 6=0.205 in "Ic", 6=0.235 
in "2c", 6=0.275 in "3c" and 6=0.320 in "4c". We also 
found approximate locations for the transitions: the "Oh" 
to "Ic" transition occurs at 6 = 0.19, the "Ic" to "2c" 
transition occurs at 6 = 0.22, the "2c" to "3c" transition 
occurs at 6 = 0.26 and the "3c" to "4c" transition occurs 
at 6 = 0.30. 

It is easy to keep the system in "cold" ("Xc") phases 
even for 6 < 6bulk- These phases are in principle 
metastable. However, in practice, for large enough iV, 
{N > 20), the "Xc" phases are very stable in a Monte 
Carlo simulation. This makes it possible to investigate 
the "Oc" to "Ic" transitions of interest also for L < 7. 
All our values for bc{L) for L < 7 have bc{L) < 6bulk- 

4. Relation to the finite temperature deconfine- 
ment transition. Suppose we studied a torus of un- 
equal sides, Lfj^. The most plausible assumption is that 
again one Z{N) will break first as 6 is increased from 
the phase where the entire Z'^{N) is preserved. Only 



now, which Z{N) breaks will no longer be arbitrary, but, 
rather, the one associated with the direction with the 
shortest L^, is selected to break first. Call this direc- 
tion /xg. The arguments from our previous paper say 
that for 6's smaller than this transition point there is 
no dependence on the parameters L^. Hence, up to the 
transition, all L^^s can be considered as infinite. But, we 
could equally well think about the with /i 7^ /xg as 
infinite, while keeping in mind that is finite. This 
puts the system in a situation considered in |l9j |. Then, 
there would be a transition as 6 is increased, even for 
finite N . This would be a finite temperature transition, 
which is first order at large N , and has a finite limit at 
N = 00, Tc [1^ . The simplest consistent assumption is 
that Ic = This is in full agreement with the viewpoint 
of the authors of l2CI. 

In reference |l8| the string tension in lattice units is 
found to behave approximately as follows: 



1 / 11 



(19) 



With a simple minded extrapolation to = 00 we obtain 
from |0| a/ctJ « 6.05. Combining this with our value of 

0.64. The most up to date 



0.26 we find 



1 



value for the infinite N value of can be found in Pi| . 

It is about 0.60. Thus, Ic — is consistent with what 
we know to date, but the evidence is not compelling. 

The special physical effects surrounding the finite tem- 
perature transition in pure YM in the planar limit were 
first discussed in [i^. More aspects have been studied 
in ji^]. Earlier numerical studies of the infinite N finite 
temperature transition in SU{N) gauge theories can be 
found in [2^ . 

5. Preserving I independence in the meson sec- 
tor. We have emphasized the volume independence of 
the pure gauge theory in the large I phase. It is natural to 
ask whether fermions moving in the gauge backgrounds 
typical to this phase also will behave as if the volume 
were infinite. This question needs to be sharpened be- 
cause we are considering here only a finite number of 
flavors, which makes the fermions "quenched" as a result 
of the large color (N) limit. The fermions simply pro- 
vide definitions for particular nonlocal gauge invariant 
observables, but do not influence the distribution of the 
gauge background. The I independence holds only for 
single traces of Wilson loops, and there is a way even for 
a very large loop to fold up into the torus. But a trace 
of the product of a fermion by an antifermion propagator 
only depends on the end points, not on a path, so there 
seems to be no way to describe a separation that would 
not fit into the torus. 

However, there is a trick which seems to allow the 
definition of fermionic observables on the finite torus 
which nevertheless describe propagation at larger dis- 
tances. This trick is at the heart of our proposed short- 
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cut to the planar limit in the meson sector. It is a sim- 
ple generalization of work in 0, [2^ . This produces a 
prescription for calculating qq properties in the large N 
limit while preserving volume independence. At the di- 
agrammatic level it is easy to understand what we have 
done jl^, but, to be sure, the arguments supporting this 
construction are far from rigorous. For this reason, we 
have undertaken an extensive test in two dimensions [2y| , 
which came out favorable. We believe that this provides 
sufficient grounds to go ahead and see what happens in 
four dimensions. 

To be concrete, let us consider the scalars M(x) = 
-^-ipxix) and M{x) = -^xi^[x). These meson fields 

are color singlets, x and y are sites on an lattice. Nor- 
mally we would expect only momenta = ^fc^ with 
= 0, 1, L — 1 to be accessible. We claim that large 
N reduction makes it possible to interpret data obtained 
on the lattice as providing predictions for momenta 
on an {NL)'^ lattice, at leading order in the expansion. 
The momenta are now written as, K -\- Q where K is as 
before, and — j^qfi with q^ = 0,1,. ..,7V— 1. We 
are after an expression for the meson-meson propagator, 
S{K + Q). We first define a shifted hnk field ul''\x) by 

C/W(2.) ^ elffi9M[/^(a;) (20) 

and denote the fermion — and X ~ X propagators 
on the lattice, in a given gauge background, {U}, by 
G{x,y, {U}). The shifted gauge field links are not in 
SU{N). Let us consider the collection of all Wilson and 
Polyakov loops made out of the ujf'' (x) variables as ele- 
mentary links. Before taking the trace all the unitary ma- 
trices are back in SU (N) and could have been obtained 
from elementary links that also are all in SU{N). This 
works because our observables are gauge invariant even 
under U{1) gauge transformations that take the links out 
of SU (N) . The dependence on the integers is remov- 
able by a Z'^{N) symmetry transformation. Thus, so long 
we treat a fermionic observables that depends on a single 
gauge field background there is no dependence on so 
long the Z^[N) symmetry is not broken. In other words, 
we have to be in the "Oc" phase. The q^ asume their 
role as momentum "gap fillers" only when we consider 
an observable that depends on different gauge field back- 
grounds. The Dirac and color indices of G{x, y, {U}) are 
suppressed. We now define the quantity 

R{x,y;{U},q) = -^Tr[Gix,y,{U^^^)Giy,x,{U<^°^)] 

(21) 

The trace in the above equation sums over color and 
Dirac indices. The key formula for S{K + Q) is given 
by: 

S{K + Q)^Y1 e"^'^""^^ {R{^, y; m, q)){u} (22) 
The averaging over {[/} restores translational invariance. 



If we are not in the "Oc" phase the procedure fails. 
One could try to impose a "Oc" phase by quenching the 
links on a 1^ lattice, but implementing the additional 
averaging puts too big a burden on the numerics and 
is likely less practicable than the procedure we propose 
here. 

6. Conclusions. Our previous work and the 
present paper make a plausible case for the following sce- 
nario: Planar QCD on a torus of side I has a nontrivial 
phase structure as a function of /. When / decreases from 
00 to any I < Ic, the system undergoes a phase transi- 
tion where the global Z'^{N) symmetry breaks sponta- 
neously. It is likely that = where is the infinite 
N limit of the finite temperature deconfinement transi- 
tions of SU{N) YM theory at finite iV's. The precise 
determination of Ic, including error estimates, and of the 
existence and locations of other continuum transitions 
needs substantially more numerical effort than invested 
to date, but is a feasible numerical project. 

A distinctive property of the I > Ic phase is the l- 
independence of arbitrary Wilson loops, which provides a 
continuum realization of lattice Eguchi-Kawai reduction. 
While the free energy does not depend on I at leading 
order in a generic, free, string theory with toroidal target 
space, among field theories, only certain gauge theories 
in the planar limit can exhibit such a property. This 
scenario, in turn, leads to the conjecture that at I > 
Ic the I dependence of meson propagators can also be 
made to disappear at leading order in the expansion 
by quenching, providing the opportunity for a numerical 
shortcut to the infinite N meson sector of SU{N) YM 
theory in infinite space-time. 
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